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Abstract 



There have been numerous studies of entanglement in spin systems. 
These have usually focussed on examining the entanglement between indi- 
vidual spins or determining whether the state of the system is completely 
separable. Here we present conditions that allow us to determine whether 
^^ I blocks of spins are entangled. We show that sometimes these conditions 

^^ . can detect entanglement better than conditions involving individual spins. 

We apply these conditions to study entanglement in spin wave states, both 
when there are only a few magnons present and also at finite temperature. 



1 Introduction 



The realization that entanglement is a resource for a number of useful tasks in 
quantum information has led to a tremendous interest in its properties, quan- 
ta I tification and in methods by which it can be produced. One area, which has 

been fruitful, is the study of entanglement in many-body systems (for a review 
see [H). 

Spin systems, in particular, have received a great deal of attention, and this 
has led to the formulation of several conditions for determining whether the 
state of a spin system is entangled. For example, in an N qubit system, with 



^. = iE4^ (1) 



fc=i 



with / = 1, 2, 3 and a^ being the PauH matrices for the k*-^ qubit, the state is 
entangled if 



(Ji)2 + ( J2>2 N 
If this inequahty is satisfied, the N-qubit state cannot be expressed as 



E 



P,f^'^ff...^pf\ (3) 

(k) . . . fli . 

where the 1 > Pj > sum to one, and Pj is a density matrix for the k qubit. 
A state that can be expressed in this form is known as completely separable. 
Further criteria for entanglement in spin systems have been developed [Sj [H [5] . 
A comprehensive study of entanglement conditions employing quantities that 
are at most quadratic in the collective spin operators has been given in [6l [7] . 

Recently two inequalities have been developed for the detection of entan- 
glement [S]. Consider a system consisting of two subsystems, which we shall 
denote by a and b. The Hilbert space for the total system is H = Ha ^S) 'Hb, 
where Ha is the Hilbert space for subsystem a and Hb is the Hilbert space for 
subsystem b. Let A be an operator on Ha and B be an operator on Hb- The 
state of the total system is entangled if either 

\{AB^)\'^ >{A^AB''B), (4) 

or 

\{AB)\^>{A^A){B^B). (5) 

Note that the condition (|4]) holds if the left-hand side is replaced by |(^t^)p. 
In the following we shall use the form (j4|) or KAt^)^ > {A^ AB^ B) interchange- 
ably. These are sufficient conditions for entanglement; if they are not satisfied 
we cannot say whether the state is entangled or not. Some of the implications 
of these conditions have been explored for the case that the systems are field 
modes [SI [9] , and in the case that one is a field mode and one is a collection of 
atoms [lOj . Here we would like to explore their implications for spin systems. 
In that case, we have a collections of spins, and our subsystems are two non- 
overlapping subsets of the total set. We describe each subset by a collective 
spin, 3a for set a and Jf, for set b. Let Ja- be the angular momentum lower- 
ing operator for set a and Jb- be the angular momentum lowering operator for 
set 6. The corresponding raising operators are Ja+ and Jb-|_, respectively. Our 
entanglement conditions become 

\{Ja-Jb+)\ > {Ja+Ja-Jb+Jb-}, 

\{Ja-Jb-)f > {Ja+Ja-){Jb+Jb-). (6) 

These inequalities differ from the ones discussed in the previous paragraph in 
that they detect entanglement between two blocks of spins, and not whether 
the state is completely separable or not. For example, if one is studying entan- 
glement in a spin chain, one may simply be interested in whether the state is 



entangled or not, in which case Eq. ^ could be of use. However, one might 
instead wish to find out whether two blocks of spins are entangled, in which 
case the above equations, with Ja and Jb being collective spin operators for the 
respective blocks, could be useful. 

As was mentioned previously, a large amount of work has been done on 
entanglement in spin systems, in particular spin chains with different types of 
interactions between adjacent spins. Typically the concurrence of two spins in 
the ground state of the system is found. This has been done for many variants 
of the Heisenberg model, both with and without an applied magnetic field. We 
will be interested in applying our conditions to a spin system at finite tempera- 
ture, in particular one with a ferromagnetic Heisenberg interaction. Wang and 
Zanardi showed that in a one-dimensional ring with periodic boundary condi- 
tions, and the spins interacting via an isotropic Heisenberg interaction, there 
is no entanglement between any two spins at any temperature [TT]. This re- 
sult depends on the SU{2) symmetry of the Hamiltonian, and does not appy if 
there is an applied magnetic field. Asoudeh and Karimipour look at the case 
in which there is an applied magnetic field and only the ground state and the 
state with one flipped spin are populated, and found the concurrence between 
two spins as a function of their separation |12] . The advantage of applying the 
entanglement conditions in the previous paragraph is that we can study the 
entanglement between blocks of spins, and not just the entanglement between 
individual spins. 

The paper is arranged as follows. We will initially study some general prop- 
erties of our entanglement conditions. We will first compare the application 
of our conditions to individual and collective spins. It will be shown that the 
conditions can detect entanglement between spins in angular momentum intel- 
ligent states. These are states for which the uncertainty relations for angular 
momentum operators are satisfied as an equality. We will then show how the 
above conditions can be strengthened by imposing local rotational invariance. 
Next, we will move on to spin systems, and use the above conditions to study 
entanglement in spin wave states, first for states containing a small number of 
spin waves, and then for states at finite temperature. 

2 Examples of states 

One of the advantages of the entanglement conditions in the Introduction is that 
they allow us to look at the entanglement between blocks of spins rather than 
between individual spins. A standard approach when studying the entanglement 
in spin- 1/2 systems is to choose two spins and calculate their concurrence. It 
is, however, quite possible that there is no entanglement between individual 
spins, but there is between blocks of spins. In that case, the method based on 
concurrence will fail. This can be illustrated by an example. 

Let us consider four qubits, i.e. spin-1/2 particles, with qubits 1 and 2 
in block a, and qubits 3 and 4 in block b. Each qubit has an orthonormal 
basis {|0), |1)}, and a raising operator ct*^"*") and a lowering operator cr^~\ where 



cr(+)|0) = |1), cr(+)|l) = 0, and o-(-) = (cr(+))t. Let us now consider the four- 
qubit state 

I*) = -l|00)i2 100)34 + ^(|01)i2|10)34 + 110)12101)34). (7) 

Tracing out qubits 2 and 4 we find the reduced density matrix for qubits 1 and 
3 

Pi3 = ^100)13(001 + i(|01)i3(01| + |10)i3(10|), (8) 

which is separable. So, if we just look at qubits 1 and 3, i.e. one qubit in each 
block, we do not see any entanglement. However, setting 

Ja- ^^['^ + O^-'^ , J6_ - 4"' + ^^\ (9) 

we find that 

{,Ja-Jb-) = ^, {Ja+Ja-){Jb+Jb-) =\. (10) 

SO that the second entanglement condition in Eq. ^ is satisfied. Therefore, 
by looking at entanglement between blocks, we see that the state is, in fact, 
entangled. 

We first will proceed to examine two more complicated types of entangled 
states in order to see whether our entanglement conditions can show that these 
states are indeed entangled. In the first case, we will apply the conditions to 
both individual and collective spins in order to see which method yields a more 
sensitive test of entanglement. 

2.1 Correlated sets of qubits 

Suppose we have 2n qubits. We will divide the qubits into two blocks of n 
qubits each, and within each block, we will consider only those states of total 
spin i = n/2. In particular, we want to examine states of the form 

I*) = Y^ Cm\3.m)a®\],m)b (11) 

m=-j 

with j — n/2, and the state with subscript a referring to the first block and the 
state with subscript b referring to the second. This is clearly an entangled state, 
and we want to see whether the entanglement conditions we have proposed will 
detect the entanglement. We will do this in two different ways. First, we will 
apply the second entanglement condition in Eq. © to the collective spin of 
each block. Next, we will choose one qubit from each block and apply the same 
condition to those two qubits. 

The calculations for the condition using the collective spins is straightfor- 
ward, and we find 

i 

(*|Ja-J6-|*) = Y. (i+"i)(i-"^+l)C-lCrn, 

T71 — — J + 1 



{^\Ja+Ja-m = {^\Jb+Jb-m= J2 |Cm|'(j+m)0--m+l). (12) 

m=-j 

Therefore, the second entangleraent condition in Eq. ([5]) becomes 



X! 0' + ™) 0' - "1 + i)<»-ic« 

m— — j + l 



> J2 Wrn\Hj+m)ij-m + l). (13) 



One possible choice of Cm is to set c,,,, = rjx^^"^, for some x > and 77 an 
appropriate normahzation constant. This gives us 



J2 ij+m){j^m+l)x^'^^+"'^-^ 



> J2 ij + m){j~m+l ) a;2 (^■+"') . 

m — —j-\-l 

(14) 
This condition is clearly satisfied when x < 1, but not satisfied for x > 1. 

Now let us see what happens if we just look at one qubit in each block. Let 
the qubit in the first block be qubit 1 and the one in the second block be qubit 
2, and we will assume that each of these qubits is the first one in its respective 
block. Let us call the spin-down state of an individual qubit |0) and the spin- up 
state |1). The basis states for each block are n-fold tensor products of spin- up 
and spin-down states for each qubit in the block. The state \j, m) of n qubits 
with j = n/2 is the symmetric linear combination of all basis states in which 

there are j + m ones and j — m zeroes. There are . I such states. The 

V J+^ J 
operator <t\ a\ , where a\ and crj are the raising and lowering operators 
for qubit 1, is just the projection onto states in which the state of the first qubit 

is 11). There are f . _, I states, and this implies that 

y J + m — 1 J 

( 2j - 1 



This implies that 



j + m 



.(+)j-)h-f,\ - sr \^ |2J 



^ar'a^'l^)^ > \cmr-^—. (16) 



m=-j 



2j 



The expression for (\1/ 1 (Tj <^2 |*) is identical. 

We now want to compute {a{j,'m—l\iSi b(j, m— l|)crJ^V2 (b, ™)a'8)|j, m)b). 
The operator a[ (J2 will pick out the basis states with ones in the first slot 



of each block. By reasoning similar to that above, we have that 

{a{j,m- 1| (g) b{j,m- l\)a[^'a^'~\\j,m)a<^ \j,m)b) 

2j-l 

j + ?71 — 1 

j + m — 1 J \ j + m 
{j + m){j - m + 1) 



This gives us that 



(*k^V(-)|vI/)= ^ 



m——j-^l 



ij + m){j -m + l) 
(2jr 



(17) 



(18) 



Finally, the second entanglement condition in Eq. ([6|) with A = a\ and B 



(-) 



0-2 becomes 



^ {j + m){j-m + l) 

m— — j + 1 



>(2j) J2 \cm\^{j+m). 



(19) 



m=~j 



Comparing Eqs. (IT^ and P^ we note that 2j' > j — to + 1 if tti > — j + 1, which 
is the entire range of the sum. This implies that at least for states of the type 
in Eq. (fTTj). the collective spin condition is stronger, that is, it will be satisfied 
by more states, than the condition for individual spins. 

2.2 Angular momentum intelligent states 

We first want to find some spin states that satisfy our entanglement conditions. 
One possibility is to find states in which the spins of the two subsystems are 
highly correlated, and states in which the uncertainty of the sum (or difference) 
of the two spins is small will satisfy this condition. 

Let us begin by looking at the uncertainty relation for the total spin, 



AiJal + Jbl)A{Ja2 + Jb2) > ^\{Ja3 + Jbs)], 



(20) 



where Jai, Ja2, and Ja'i are the components of Ja, and Jf,i, Jm, and Jbs are the 
components of Jf,. We would like to find the states which satisfy this relation 
as an equality. These states were first found in [TB] , and here we will follow the 
treatment given in |14j . These satisfy the eigenvalue equation 



[(Jla + ^16)+iA(J2a + J2b)]|* 

where A is real. This equation implies that 

m{Jal + Jfcl)|*) = Re(/3), (*|(Ja2 + MW 



/?!*), (21) 

(l/A)fm(/3), (22) 



and 

[AiJal + Jtl)]^ ^ ^{Ja3 + Jbs) [AiJa2+Jb2)?^^{Ja3 + Jb3). (23) 

From these equations, we see that when A is smah, Jai and —Jbi are highly 
correlated, and when it is large, Ja2 and —Jb2 are highly correlated. These 
states are spin analogs of two-mode squeezed state for light. 
In order to solve Eq. (|2T|) . we first define a state 

I*') =e'''(^"i+^")|*), (24) 

and insert the resulting expression for |$) into Eq. ((2T|) to give an equation for 
I*') 

{{Jai + Jbi) + i\[{Ja2 + Jb2) COS - (J„3 + ^w) sin0]}|*') = /3|*')- (25) 

Now choose A = —1/ cos^, and 9 to be in the range tt > > tt/2, which implies 
that A > 1, and 



[{Ja- + Jb-) - iv/A2 - 1(J„3 + JhsW) = /?!*')• (26) 

We now expand I^I'') as 

i 
1*')= J2 Cnra\n,m), (27) 

where we have set \n,va) = |j, n) ® |j, ?ti). If we assume, for simplicity, that 
Cn,m = 0, unless n — m, our equation for |^') reduces to the recurrence relation 



f3 + 2mi\/A2 - 1 

Cm+l,m+l = T— , -|W . ■rCm,m, m < j, (28) 

[J +m + l)(j - to) 



[/? + 2ij v/A2 - 1] Q,, - 0, TO = J. (29) 

From the second equation, we see that either /3 = —2ij\/X^ — 1, or C^j- = 0. 
If Cj,j = 0, then it must be the case that (3 = —2moi\/X'^ — 1 for some toq. So, 

Cm^^ = i-2^V^-ly-- ^'^^^^'Zi-. C^.-. (30) 

(mo-TO)!(j +to)!(2j)! 

for TO < TO-o and Cm™ = for m > Trig. After grouping the m- independent 
constants ^^^§7]! ^-j,-j i^^^o Cj,m(,(A), we have 

|*(j,mo,A)) - Q,™„(A)e-^''(^-+-'") 

Y, {-2jV>F~iy+^- — ^L_^^i--^|TO,TO). (31) 

^-^ . ymf) — my\] + ray. 



We want to see if there is some range of parameters for which \{Ja-Jb-)\^ > 
{Ja+Ja-){Jb+Jb-), but for thesc states (Ja+Ja-) = (Jb+Jb-), so we just need 
to show that \{Ja-Jb-)\ > {Ja+Ja-}- We find that 



\Ja+Ja- 



|C,,™o(A)P Y. [4(A^-l)F+™[ 



(j - m)\ 



{nio — my.{j + my. 



(32) 



and 



(Ja-Jb-) = \C,..n„iX)\' Y. [4(A'-1): 



j+rn r 



{j -my. 



2iVA2 + 1 
A 



(rng — m) 



{niQ — m)!(j + m,y. 

(A^-l)m^ " 
A2 



(33) 



Consider the simple case in which j — l^rriQ — — l,and tti = —1, so that the smu 
has only one term. The entanglement condition becomes 



A2-1 



A^ 



and the state is entangled if A > 3. 



> 



(A + 1)^ 
2A2 



(34) 



3 Local rotational invariance 



Entanglement is not affected by local unitary transformations, and so, ideally, 
we would like our entanglement conditions to be invariant under local unitaries 
as well. It is not always possible to accomplish this, but we can sometimes obtain 
invariance under a subgroup of the group of local unitary transformations. In 
Ref. [TU] in which entanglement between field modes was considered, it was 
possible to find entanglement conditions that are invariant under local Gaussian 
transformations of the field modes. These new conditions were stronger than the 
original ones, that is they detect entanglement in a larger set of states. Thus, 
making the conditions invariant under a subset of local unitary transformations 
strengthens them. 

For the entanglement conditions we are considering in this paper, the obvious 
group of local unitaries consists of local rotations. Under the action of the 
rotation _R(a, /3, 7) 



^,,aJi^-t0j2^-tfJ3 ^Q have that 



1 i 

R^^Jj^R — [-(cosa + cos/3) + - sinasin/3]e*''J4_ 

1 i 

+ [-(cos/3 — cosa) H — sinasin/3]e *'''J_ 

+ [1 — i sin a sin /3] e*^ J3 . 



(35) 



Now suppose we start with the entanglement condition Eq. Q, and we want to 
find from it a condition that is invariant under local rotations of the a system 
(finding a condition that is invariant under local rotations of both a and b 
subsystems is possible, but it results in a 9 x 9 matrix, which is rather unwieldy). 
We note that what the local rotation on subsystem a does is to send both Ja+ 
and Ja- into linear combinations of Ja+, J a--, and Ja^- This suggests that we 
set A = c\Ja- + C2 Ja+ + C3 Ja3 and B ^ Jb_ in Eq. ^. The entanglement 
condition can then be written in the form 

{cl cl c^ ) Af j C2 > 0, (36) 

where M is a 3 x 3 matrix, whose elements are linear combinations of expectation 
values of products of angular momentum operators. In particular, 

Mil — \{Ja-Jb+)\ — {Ja+Ja-Jb+Jb-), 

M12 — {Ja-Jb+)*{Ja+Jb+) — {Ja+'h+Jb-), 

Mi^ = {Ja-Jb+)*{Ja3Jb+) — {Ja+JasJb+Jb-}, 

M21 = {Ja+Jb+)*{Ja-Jb+) — {Ja-Jb+Jb-), 

M22 = \{Ja+Jb+)\ — {Ja-Ja+Jb+Jb-), 

-^^23 — {Ja+Jb+Y {JazJb+) — {Ja-JazJb+Jb-), 

M^i = {Ja3Jb+}* {Ja-Jb+) — {JaiJa- Jb+Jb-) , 

M32 = {Ja3Jb+)* {Ja+Jb+) ~ (JasJa+Jb+Jb-), 

M33 = |(Ja3Jb+)P - {Jl^Jb+Jb-). (37) 

If we change the state by a local rotation of system a, the effect on Eq. (|36p 
is only to change the values of Ci, C2 and C3. This follows from the fact that 
when A is conjugated by the rotation Ra, the form of the operator stays the 
same, that is, it is a linear combination of Ja+, Ja-, and Jas-, but the coefficients 
multiplying the operators change. If the matrix M has a positive eigenvalue, 
then we can find values of ci , C2 and C3 so that the above condition is satisfied, 
simply by choosing them to be the components of the vector corresponding to 
the positive eigenvalue. Therefore, our new entanglement condition becomes 
that M has a positive eigenvalue, and this condition is invariant under local 
rotations on system a. 

Let us show that this new condition is stronger than our original condition. 
If the state we are considering is 

l*> = 4s(l - J-- -J + 1> + I - J + 1> -j)). (38) 



then 



V2^ 



M = -2p . (39) 



3' 











-2p 











-3' 



Noting that j^ is positive, we see that the state is entangled. Because this 
condition is invariant under rotations of system a, it would also show that the 
state Ra (g) /blvP) is entangled. 

Now, let us see what happens if we apply our original condition to the state 
Ra(^ hl"^)- We begin by finding 

2 ,-2 

{^\R-\Ja+RaJb-\'i') = ^[icosa + COS 13)^ +sm^ a sin^ P], (40) 

and 

o2 

i^lR-^Ja+Ja-RaJb+Jb-l'i') = ^ [(cos ^ - COS a)^ + siu^ a siu^ /3] 

+j'''(l + sin2asin2/3). (41) 

Therefore, the state is entangled according to the old condition if 

cos a cos /3 > j(l + sin^ asin^^) + -(cos a cos /3 — 1)^. (42) 

This condition can be satisfied for only a limited range of a and /3 if j is small, 
and it cannot be satisfied at all if j > 1, which actually allows only for j = 1/2. 
Therefore, our new condition, which is invariant under rotations of system a, is 
considerably more powerful in that it detects entanglement in a much larger set 
of states. 



4 Spin waves 

The low- lying energy states of a system of spins coupled by exchange interactions 
are wavelike, as shown originally by Bloch for ferromagnets. The waves are 
called spin waves, and they correspond to excitations of definite energy called 
magnons. We will study the entanglement between spins, and blocks of spins for 
magnon states in a ferromagnet. We will first examine entanglement in states 
containing a small number of magnons, and then go on to study the case of a 
ferromagnet at low, but finite, temperature. 

4.1 Small number of magnons 

The Hamiltonian describing spins on a lattice interacting via a nearest-neighbor 
exchange interaction and an externally applied magnetic field is |151 116] 



where the vectors (5 connect the spin at site j with its nearest neighbors on 
a bravais lattice, J is the exchange integral, which is assumed to be positive, 
Mo — (3/2)/is is the magnetic moment of the atoms, and Sj is the spin angular 
momentum operator of the atom at j. Hq is the intensity of a static magnetic 



10 



field directed along the z axis, and we will take the limit as Hq — > 0+ to make 
the magnetic moments line up along the positive z axis when the system is in 
the ground state |ri). The z component of the total spin, S^ = X), ^j^ i^ ^ 
constant of the motion, and the ground state of the system simply has all of the 
spins pointing in the +z direction. 

For the case of a small number of spin waves, let us consider a line of iV 
spins with periodic boundary conditions (the spin at iV + 1 is the same as the 
spin at 1). If the atoms have a spin of 1/2, then the state containing a single 
magnon is a linear combination of states with one spin flipped 



1 ^ 
l*) = ^^E^'''"'^i"^l")' (44) 



AT 

where a is the spacing between spins, and k = 2'Kn/{Na), where n is an integer 
in the range —{N/2) < n < {N/2). The operator a^~' is the spin lowering 
operator for the spin at site j, that is it maps the spin up state at site j to the 
spin down state at the same site. 

Now let us examine the entanglement of this state using Eq. ^ . Let 

m 

L-\-m 

B = S2+= J2 ^r^' (45) 

j=L + l 

where m is a number such that 2m < N . This will allow us to see if there is 
entanglement between two blocks of spins each of size m and distanced from 
each other by {L — m) spins. Our state is entangled if 

\{Si+S2-)f > {S1-S1+S2-S2+). (46) 

For the single magnon state above, the right-hand side is zero, so as long as the 
left-hand side is non-zero, we can say that the blocks of spins are entangled. In 
fact, we find that 

-. m L-\-m 

(•^1+^2-) = ^E E e^'"^''-''^- (47) 

jl = lJ2 = L + l 

If the size of the blocks is small compared to the wavelength of the spin wave, 
the term in the sum will all have approximately the same phase, and will add 
coherently. This would show that the blocks of spins are entangled for this state. 
If we want to look at more than one magnon, more sophisticated techniques 
are required. We will make use of the Holstein-Primakoff transformation, which 
expresses the spin operators in terms of boson creation and annihilation op- 
erators, and allows us to approximately diagonalize the Hamiltonian. The 



11 



Holstein-Primakoff transformation of the spin operator Sj to boson creation 
and annihilation operators a-, Oj is given by 

5j+ = Sj^ + iS^y = {2S - a]aj)i/2„j, 
Sj- = Sj:r - iSjy = a] (25' - ajaj)^/^ 
5j, = 5-a]aj, (48) 

where 

[«j,4]-'5j,,. (49) 

If we consider only situations in which the number of flipped spins is small 
compared to the total number of spins, we can expand the square roots and 
keep only the first terms in the expansion. In addition we make a transformation 
from the spin operators, a- and Oj, to the niagnon variables, 6j^and 6k, defined 

by 

6k = 7V-i/2^e*'^""Jaj, (50) 

j 

where Tj is the position of spin j. The magnon operators satisfy boson commu- 
tation relation: 

[&k,&L] = '5k,k', [6k,6k']=0. (51) 

When the number of flipped spins is much less than N, the Hamiltonian is 
diagonal in the magnon operators, 

i? = ^c^kfo],6k, (52) 

k 

where 

wk = 2 Jz5(l - 7k) + 2^ioHo, (53) 

and 

7k = ^E^"-'- (54) 

As was mentioned before, we will work in the limit Hq -^ O"*", so that in the 
ground state the spins are lined up along the z axis. In these equations a center 
of symmetry is assumed so that 7k = 7-k, and z is the number of nearest 
neighbors each spin has. 

Now we are in a position to consider the two-magnon state. We shall again 
consider the one-dimensional case, that is TV spins in a line. We want to study 
the entanglement of the state 

I*) = blbljO) - ^E^"^'""e"''^'"«i«i|0) (55) 

for fci =/= k2. We shall examine the entanglement between two blocks consisting 
of m spins each, one beginning at spin 1 and the other beginning at spin L, so 



12 



that the blocks are separated by L — ttt, spms. Therefore, we choose 

m L-\-m 

A^V2S^aj, B^V2S Y^ a^, (56) 

j=l ]=L+1 

in Eq. (g]). We find 

492 
(A^AB^B) = — r{2xy + 2:Eycos [La{ki ~ fcz)]}, (57) 

where x = [cos(fci ma) + l]/[cosfcia + 1], and y = [cos{k2ma) + l]/[cosfc2a + 1], 
and 

\{A^B)\^ = -^{x^ + 2xycos [La{ki - A^a)] + y^}. (58) 

Therefore, the state is entangled if 

{x - yf > 0, (59) 

which is true as long as x ^ y. One situation where a; = y = 1 is when the block 
size is one m = 1, implying that the condition Q does not detect entanglement 
between individual spins in the two-magnon state. Recall that kia = TT2ni/N 
and ^20 = TT2n2/N. If the block size m is such that m2ni/N — 2li ^ \ and 
m2n2/N = 2I2 + 1 where li and I2 are integers, hence x — y — and no 
entanglement is found according to the inequality ([5^. The condition ([5^ 
indicates that in the ideal zero-temperature two-magnon state, entanglement is 
found regardless of how far the two blocks are separated. This no longer occurs 
in the more realistic non-zero temperature state we are going to investigate 
below. 

4.2 Finite temperature 

Now that we have seen that the entanglement condition, Eq. ([4]), is useful in 
detecting entanglement in states consisting of a few magnons, let us see whether 
it can also detect entanglement in a system of ferromagnetically interacting spins 
at a finite temperature, T. The density matrix for the system is now given by 

P=^e-P", (60) 

where (3 = l/(fcsT) and kg is Boltzmann's constant. The partition function of 
the system, Z is given by 

Z ^ Trie-n -HE e-^"^"^ - n JT^^y (61) 

and rik is the number of magnons with wave vector k. 

We first look for entanglement between two individual spins having radius 
vectors ri and r2 by employing the inequality in Eq. (j4]) with 

A ^ Sj,+ = V2Sai, B = Sj^,+ = V2Sa2. (62) 
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We then have that 






N 

ki k2 



Using the relationship X]riLo("' ^" ^)x" = 1/(1 ^ 2;)^i one obtains 

^^(^•^^^l^)-!^!^' (64) 

so that 

OC. . „-iki-(ri-r2) 

(ABt) = i^yf . (65) 

ki 

This gives us the left-hand side of our inequality, and we now need to find the 
right-hand side. Using 

("-ki , nk2, nk3, nk4, . . • l&k^&ka^kj^kil'^ki, ?ik2, "ka, »^k4, • ■ •> = 

(5ki,k2^k3,k4?T-kin-k3 +(5ki,k4 5k2,k3»T-ki(?lk2 + 1), (66) 

we obtain for the right hand side of Eq. ^ 



2S\ 

'n) 



ki,k2 



g-^("fci+"fc2) gi(ki-k2)-(ri-r2)g-/3Ljfcj 



(1 - e"'3"^i)(l - e-^'^"'^-) (1 - e-'^'^'=i)(l - e"'^'^'''^) 



(67) 



Equation ([57]) shows that (Atyl_Bt_B) can be separated into two parts: the first 
one represents the self correlation of the particles and is distance independent, 
while the second represents interparticle correlations and depends on the dis- 
tance between the particles. 

Let us now specialize to a cubic lattice with lattice constant a and z = 6. If 
|ka I < 1, then 

1 — 7k = 1 — -(cosfca^a -I- coskytt + cosfc^a) ~ -k^a^, (68) 

o o 

and the magnon energy can be expressed as ujk = Dk^, where D = 2JSa^. To 
tackle the sums over kj we note that 

- - < fc,- < - (69) 

a a 

and approximate the cube by a sphere, so that kj < ^/3n/a. We replace the sums 
by integrals in a spherical coordinate system, and our entanglement inequality, 
Eq. (JU becomes, upon using Eqs. (|65)) and (|67l) and carrying out the angular 
integrations, 

Q = I^-ili+hh)>0, (70) 
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where |(ABt)|2 = /f, {A^ AB'^ B) = /| + hh, 



Jo \yo a J 1-e y 

^' = 1 ^^1^^^' (^2) 

i.-hyf(^^^V3)f^, (73) 

Jo \yo a J 1-e y 

and f{x) is the famihar function 

fix) = ^^. (74) 

a; 

Here Ar = |ri — r2| is the interatomic distance, j/o = yf^D^fZ-n ja and the 
dimensionless integration variable y is related to the wave vector component k 
by y = \fpT3k. Due to the presence of the exponentially decaying factor e^^ in 
the numerators of the integrands in I2 and /a, small values oiy, y ^ 1, contribute 
most to these integrals. In the case of I3, the fact that ?/^e~*''/(l — e^y ) is a 
decreasing function, causes that integral to be positive. 

As T increases, the upper limit of the integrals, yo, which is proportional to 
-^, tends to zero, and, as a result, e~^ — )• 1 and /a — >■ /i. Hence the inequality 

([70)) becomes if — (/| + if) > 0, which cannot be fulfilled. This means that our 
condition does not show the existence of entanglement in the high temperature 
limit, which is consistent with what we expect on physical grounds, i.e. that 
there is no entanglement at high temperature. As the temperature decreases, the 
upper integral limit yo increases. The integrand in /i is an oscillating function 
of y with a varying sign and an increasing magnitude, and the sign and value of 
Ji are determined mostly by the contribution near yo- For short distances and 
low temperatures, the absolute value of I2 is typically much larger than those 
of Ji and I3, which makes it the leading factor in deciding the sign of Q. 

In Fig. [1] we give a representative example of the distance dependence of 
Q. Positive values of Q indicate entanglement. It can be seen from Eq. (|7ip 
that Ii is an oscillating function of the interatomic distance Ar, with a damping 
envelop. This shows up in the behavior of Q: If we allow for continuous values 
of Ar/a we will see the damped oscillations more clearly. For short interatomic 
distances, entanglement is clearly observed. Q turns negative for the first time at 
Ar/a = 13. However, it can again become positive, meaning the reappearance of 
detectable entanglement at much larger distances before becoming permanently 
negative. In Fig. [T] the temperature is fixed. For lower temperatures, the 
shortest distance at which Q is found to be negative and the overall range over 
which Q is found to be positive increases. 

The temperature dependence of Q is illustrated in Fig. [5] for different values 
of the inter-particle distance. It can be seen that as the temperature increases, 
Q monotonically decreases, and the shorter the interparticle distance, the later 
Q crosses into the negative range. In other words, as we would expect on 
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Figure 1: The quantity Q as a function of the interatomic distance, scaled by 
the lattice constant. A positive Q indicates entanglement. Some large (positive) 
values of Q are beyond the scope of the figure. The temperature is fixed at 
^2JS/{kBT) = 7. 




1 2 

T (units of 2JS/kB) 



Figure 2: The quantity Q as a function of the temperature for different inter- 
particle distances Ar/o = (a) 1 (solid line), (b) 3 (dashed line), (c) 10 (dotted 
line), (d) 20 (dot-dashed line). 
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physical grounds, lower temperatures and shorter inter-particle distances are 
more favorable for entanglement generation, and this parameter region is where 
our condition shows the presence of entanglement. If we assume some typical 
parameters for ferromagnets [T5] D ~ 0.5 x 10~^^ erg cm^ and a ~ iA, using 
A;b = 1-38 X 10~^^ erg K"'^, the temperatures at which Q turns negative, which 
is where entanglement is no longer detected, are 420K, 150K, 20K, and 8K for 
Ar/a = 1 (solid line), 3 (dashed line), 10 (dotted hne), and 20 (dot-dashed 
line), respectively. Thus when the atoms are closer located, entanglement can 
be detected at higher temperatures. We are, of course, assuming that these 
temperatures are still considerably below the critical temperature, so that the 
spin- wave description remains valid. 

Let us now proceed to use the condition (|4]) to investigate entanglement 
between blocks of m spins each, one beginning at spin 1 and the other beginning 
at spin L. With A and B chosen as in Eqs. (1551) . calculations similar to the 
derivation of Eqs. ([7n)) -(rM l) show that the entanglement condition now takes on 



the form 



m L+jn \ 2 f' r mm 



3ii' 
m L-\-m \ / 7n L-\-7n \ \ 

E E ^-.)(E E ^3.) >0, (75) 

where luj and hij are given by the respective Eqs. ([7T|) and ([73]) with Ar = 
|ri — r2| being replaced by Ar^ = jr^ — Tj\. Again {A"^ AB'^ B) (the term in 
the curly brackets) consists of two parts, the first representing the correlations 
between spins within a block and the second representing the correlations be- 
tween blocks. Equation (j75l) is general in that the atoms can be arranged in an 
arbitrary manner in space. The only assumption used is that the two blocks do 
not overlap. As in the case of individual spins, in the high temperature limit 
2/0 — >■ 0, I'iij — > Iiij, indicating explicitly that the inequality ([75)1 cannot be 
satisfied. At low temperatures, whether Q is positive or negative depends on 
the details of the terms in the sums over luj . 

In Fig. [3] we plotted Q as a function of the block size to, to, being the 
number of spins contained in each block. It is assumed that each block consists 
of neighboring spins located along a straight line, one beginning at spin 1 and 
the other beginning at spin L so that the blocks are separated hy L — ni spins. 
For the parameters used in Fig. |21 the case of individual spins to = 1 exhibits no 
entanglement (cf. Fig. [1]). As the size of the blocks to increases (Fig. |3l inset), Q 
acquires positive values indicating a presence of inter-block entanglement. The 
change is not monotonic, however. As to increases further, the entanglement 
detected by our condition can disappear and reappear, being particularly strong 
for TO = 11, 12, 13. The sign of Q obviously depends on whether the luj add 
constructively or destructively. An examination of inter-block entanglement 
may thus offer much richer physics than simply a study of the entanglement 
between individual spins. 
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Figure 3: The quantity Q as a function of the block size m - the number of 
spins in each block for a fixed dimensionless temperature ^/2JS/(k^^T) = 7 and 
for L — 13. The inset zooms in the part of the plot for small m. 

5 Conclusions 

We have presented two entanglement conditions for spin systems that allow us to 
study the entanglement between blocks of spins. Most tests for entanglement in 
spin systems test for either complete separability or for entanglement between 
individual spins, and the results in this paper complement those. We have 
shown that in some cases the conditions involving blocks of spins can detect 
entanglement when tests of individual spins cannot. It was shown that our 
entanglement conditions can detect entanglement in intelligent spin states and 
in states of a spin chain containing a small number of spin waves. This latter 
result was then extended to show that entanglement in spin waves at finite 
temperature can also be detected. 
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